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The influence of the homogeneous magnetic field on a single mobile hole in a magnetic insulator, 
as represented by the two-dimensional t-J model, is investigated by considering the coupling of the 
field to the orbital current. The energy of the J = system is analyzed via the high-temperature 
expansion and the small system diagonalization. The susceptibility is shown to be diamagnetic and 
diverging at low temperatures T. In contrast, in the antiferromagnetic J > case small systems 
generically reveal a tendency towards a paramagnetic response in larger fields at low T. By employing 
at T = the cumulant expansion we study the ground state in arbitrary B, showing a behavior 
very sensitive to the character of the quasiparticle dispersion and the magnetic-field strength. At 
low B the perturbation and small-systems results are consistent with a pronounced diamagnetic 
susceptibility at T ^ 0, but indicate on a suppressed contribution at intermediate T ~ J. 



PACS numbers; 71.27.-f-a, 75.20.-g, 74.72.-h 
I. INTRODUCTION 

In a system of correlated electrons the external mag- 
netic field can induce several interesting effects. Theoret- 
ically the least understood are those phenomena, where 
the magnetic field couples to the motion of charge car- 
riers. In recent years it has been realized that the un- 
derstanding of the anomalous temperature- and doping- 
dependence of the Hall effect is one of the most chal- 
lenging questions in connection with the normal state of 
cuprates, being representatives of two-dimensional (2D) 
doped magnetic insulators. Here even the sign of the ef- 
fect is theoretically controversal |^,||. The diamagnetic 
contribution to the d.c. susceptibility has been much less 
investigated M, although it is closely related to the Hall 
conductivity B- It emerges from the orbital motion of 
mobile carriers. For noninteracting electrons the contri- 
bution corresponds to the Landau diamagnetism, which 
is largely temperature independent. In analogy to the 
Hall effect and other anomalous properties of the nor- 
mal state in cuprates, one could expect anomalies also 
in the diamagnetic contribution. So far, however, both 
experimental and theoretical answers are lacking. 

Magnetic field dependence of the eigenstates of tight- 
binding electrons is nontrivial even in the absence of any 
electron correlations 1^, in particular when the depen- 
dence of the ground state on the field strength B and elec- 
tron density is investigated 0,1). There have been only 
few analogous studies of correlated systems. Recently, 
the ground state of a single hole in the 2D t-J model in 
the presence of magnetic field has been studied. The 
main message is that for finite (but not very small) B the 
energy is reduced by an amount proportional to parame- 
ter t, and the result in such a doped insulator was inter- 
preted in terms of the composite nature of quasiparticles 
(QP) [^,0. Another evident observation is, however, 
the difficulty to extract a reasonable result from studies 



of small systems. 

The aim of this paper is to elaborate on the problem 
of a single hole in a magnetic insulator in the presence of 
a homogeneous magnetic field. We study the planar t-J 
model as a prototype model for strongly correlated 
electrons and electronic properties of cuprates, 



+ R.C.) + JJ2^^■^J. (1) 



where c. 



are fermionic operators, projecting out 



states with the double occupancy. We consider the sys- 
tem in a homogeneous field B, perpendicular to the 
plane, and use for convenience the Landau gauge, where 



9y = ■^d,j--A(r,), 



B{0,x,0), 
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with d, 



The relevant parameter for the 



strength of B is the dimensionless flux per plaquette 
a — iTiBa^/ where 0o = h/e is the unit quantum 
flux, and the relevant regime is — tt < a < tt. Further on 
we set the lattice spacing oq = 1, as well as = fc^ = 1. 

In the following we restrict our study of the model 
Eq. (0) to the case of a single hole doped into a magnetic 
insulator. The idea is that results for a single hole (spin 
polaron) remain relevant for the regime of finite, but low, 
hole concentration Ch 1. Here a semiconductor-like 
picture implies, so that most measurable quantities, as- 
suming the independence of spin polarons, should simply 
scale with Ch- E.g., the diamagnetic susceptibility should 
behave as x k c/j. 

The ground state of the spin polaron at _B = has been 
studied extensively both by analytical and numerical ap- 
proaches, and can be considered as one of few rather 
settled problems within the theory of correlated systems. 
Still, here at least two substantially different regimes have 
to be distinguished. 

At finite J > (as relevant for cuprates with J /t ^ 0.3) 
the ground state of a hole in an antiferromagnetic (AFM) 
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spin background has the property of a quasiparticle (QP) 
with 5=1/2 and a weU defined dispersion eo(k). Con- 
sistent results have been obtained for eo(k) using the 
self-consistent Born approximation (SCBA) ||l^,|l^, per- 
turbation expansion , numerical approaches including 
both the exact diagonalization of small systems and the 
quantum Monte Carlo method [Q. Calculations repro- 
duce a minimum at k* — {±1: 12,^7: /2), which is very 
anisotropic, i.e. /i ~ mi_/m\\ 5 for J/t ^ 0.3. This 
indicates a very weak dispersion along the AFM zone 
boundary, connecting k = k* with k = k** = (tt, 0), 
(0,7r). Studying small systems Q the latter dispersion 
is not easy to reproduce correctly. E.g., on a frequently 
studied system of 4 x 4 sites, states with k* and k** are 
degenerate, so pronounced finite size effects are expected. 
Since a small B just probes the effective mass of the QP, 
it is not surprising that results obtained on small lattices 
are not reliable or can be even misleading On the 
other hand, ARPES measurements on undoped cuprates 
[^6| show a more isotropic minimum around k* . The ex- 
planation seems to be beyond the simple t-J model, and 
the additional effect is attributed to the next-nearest- 
neighbor hopping (t') term ]l7| |. 

The behavior at J = is quite different. As shown 
by Nagaoka |l^] , the ground state is ferromagnetic (FM) 
with S = 5'inax and momentum k = 0, where the QP is a 
simple hole in the filled band of polarized electrons with 
an unrenormalized band mass. Nevertheless, close to this 
simple QP branch there is a large density of complicated 
excited states, where the hole motion is predominantly 
incoherent ||T^ . Therefore it is expected that even moder- 
ate temperature T > should have a considerable effect. 

The paper is organized as follows. Section II is devoted 
to the study of a single hole at J = and arbitrary B. 
Results are obtained via the high-T expansion and the 
Lanczos diagonalization technique for small systems at 
T = and for T > as well. In Sec. Ill we consider the 
AFM case with J > 0. Here the analysis of small systems 
at finite B is employed together with the study of the 
ground state using the cumulant expansion int/ J, which 
has proven to be very informative for B — |lj]. In 
the last section. Sec. IV, our results are summarized and 
a brief discussion on the magnitude of the susceptibility 
and its relation to the Hall constant is presented. 



II. J = CASE 



F = -TlnZ = -TlnTri 
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is within the high temperature expansion expressed in 
terms of moments /i„ and cumulants A„ 



\nZ = hiTrl + ln 



n=l 



1 + E 

InTrl + y ^A, 

r?,l 
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where f3^t/T, H = H/t and 

^„ = (_l)"Tr7?'7Trl. 
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Moments fin can be expressed as a sum over (^"2) 
closed graphs (paths). Counting different spin configura- 
tions for B — which remain unchanged jl^ after the 
performed path, each graph contributes a weight 2^~^^^ . 
Here / is the number of cycles in the spin permutation 
resulting from a hole traversing the graph, and r is the 
number of different sites in the graph. For _B > the 
only change comes from the contribution of the enclosed 
magnetic flux, so that the weight becomes 



(6) 



where m is the area of the graph in units of a§ . 

Here it is helpful to choose a 45° rotated coordinate 
system so that a 2D graph decouples into a direct prod- 
uct of two ID graphs. In this way it is straightforward 
to generate nonequivalent graphs numerically. We were 
able to evaluate /i„ and A„ up to the order n = 18. In 
Table | lowest cumulants (n < 6) A„ = ^rnn cos ma 
are presented for illustration, while higher cumulants are 
available upon request. 



A. High-T expansion 
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To study a single hole, as described by the model 
Eq.(|l]) with J — Q and _B > 0, we first use the stan- 
dard high-T expansion (HTE). Its application is in this 
case simple, since the only expansion parameter is t/T, 
while B remains arbitrary. The free energy 
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Prom Eqs.(^,||j|) the series for the orbital susceptibihty 
(per one hole) can be generated, 



X = fJ-o 



Xo 12'^ ^ 120'' 



(7) 



B=0 



2087_^, 



16128' 



8161 ^9 
53760^ 



(8) 



where xo = fJ-oe'^o.ot/fi^- 

There is no unique procedure for the extrapolation of 
the power series Eqs.(^,||) to low T. For the present prob- 
lem the most reasonable approach seems to be via the 
density of states p{e) and their moments |l9[| , 



Z ^ J p{e) e-^'de, 



Mr, 



£'>(e)de. 
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The density of states can be expanded in terms of Leg- 
endre polynomials, 



with coefficients 

Cn 



r(t + i)r(t + i) 
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for n > £ and even £ + m, while Cni = otherwise. 

The density of states is used to extrapolate both F and 
X to low T. After solving the linear equations ( p^ ) for 
Bi, we can calculate the susceptibility 
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where the same set of equations as in 

(|-[l]) holds for 

derivatives with respect to a as well. 



B. Small system diagonalization 

In analysis of the t-J model the Lanczos technique for 
the exact diagonalization of small systems has been al- 
ready extensively employed predominantly in the 
investigation of the ground state and their static and 
dynamic properties. Recently a method combining 
the Lanczos procedure and the random sampling has 
been introduced |2^ which allows calculation of finite- 
temperature properties in small correlated systems. The 
method has been used in the study of various response 



function within the t-J model at T > Ig^l- The appli- 
cation is particularly simple for static quantities, which 
can be expressed as expectation values of conserved quan- 
tities l^lj]. The calculational effort is comparable to the 
ground state evaluation. In particular the average energy 



K M-1 
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N. 
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K M-1 
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should be evaluated in this way, where |V'm);^m ^6- 
spectively approximate eigenfunctions and energies ob- 
tained by the diagonalization within the orthonormal set, 
generated from the initial functions |n) in M Lanczos 
steps. The K initial functions \n) are chosen at random, 
while Nst is the dimension of the complete basis. Note 
that it is enough to choose M, K ^ Nst ■ For more de- 
tailed explanations we refer to Refs. pO| , ^ . 

The introduction of finite B > in the model, Eq. (P , 
reduces the translational symmetry and thus for a given 
system size increases the required minimal basis set. We 
are at present able to consider the problem of a single 
mobile hole in canonical 2D systems with N = 16, 18, 20 
sites and periodic boundary conditions (p.b.c). 

It is nontrivial to incorporate phases due to a homoge- 
neous -B, being at the same time compatible with p.b.c. It 
is well known that this can be accomplished only for 
quantized magnetic fields B = ttiBq, where Bq = ipo/N 
is the smallest field corresponding to the unit quantum 
flux per system. To incorporate such B in small sys- 
tems (tilted squares), the following procedure is used: 
a) phases Oij for all hops inside squares are left as given 
within the particular Landau gauge, b) phases attributed 
to hops across the square boundaries are subject to the 
condition that the magnetic flux in each plaquette remain 
the same, i.e. B' = B mod((/!'o) (up to the addition of an 
unit flux per plaquette). These boundary requirements 
lead to a set of linear equations which have solutions only 
for B = ttiBq. Only for an untilted square lattice with 
N = L X L sites the phases can be expressed in a simple 
form as 



^(icc,iy)(icc + l:iB) 
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ix < L, 
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C. Results 



(14) 



Let us first discuss the polaron internal energy {e{a)) 
as obtained from the high-T expansion 
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From cumulants \nm in Table I we see that the a de- 
pendence of {e) first enters within the order [3^ . Such a 
term originates from a hole hopping around a loop, con- 
tributing to when all spins are equally polarized, 
in analogy to the processes contributing to the Hall con- 
stant [Q. As a result, a dependence of (e) is vanishing 
fast (oc [t/T f ] for T > t. 
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FIG. 1. Polaron energy (e) (in units of t) for J = at two 
different temperatures, as a function of dimensionless mag- 
netic field a. Full line represents high-T expansion results, 
while open and full dots correspond to results on systems 
with A'^ = 16, 18 and 20 sites, respectively. 

In Fig. |l| we present (e(a)) for lower T, i.e. T — t and 
T = 0.2 t. The results of the high-T extrapolation and 
finite size calculations performed for N — 16, 18 and 20 
are compared. For most points the agreement between 
both methods is quite satisfactory. On the other hand, 
there are clearly visible anomalies at a = wk /2 for N = 
16 and less pronounced at a = for = 18, 20. It 
is straightforward to explain these discrepancies within 
the high-T expansion. In small systems there are some 
additional processes due to p.b.c, which lead to changes 
of (e(a)) relative to an infinite system. E.g., in the N — 
4x4 system contributions from graphs in lowest field- 



dependent order /3^, representing four consecutive hops 
in the x- or y-direction cancel, except for a = n7r/2, 
changing the cumulant 



AA4 = 1, 



(16) 



as seen in Fig. |[ 

It should be observed that for — 16 this correction 
to (e) is within the leading order of [3'^, while analogous 
corrections in larger systems, e.g. for N = 18, emerge 
only in higher orders. This confirms that on small sys- 
tems the calculation of i?-induced diamagnetic currents 
is more delicate than the evaluation of most of the static 
polaron properties. Nevertheless, at J finite size ef- 
fects are rather well under control, at least in comparison 
to the AFM case J > presented in Section III. 

Consistency of high-T expansion and small system 
(A^ > 16) results allows reliable extrapolation of the sus- 
ceptibility X to quite low T ~ 0.1 t, using the procedure 
via the density of states p(e), Eqs.(p[jl^). The result is 
presented in Fig. ^ and as expected x is diamagnetic. 
While for T ^ t the susceptibility x is proportional to 
/3^, the variation is less steep for 0.1 < < T < t, where 
the variation is closer to x oc /?'' with 77 < 1. 
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FIG. 2. Susceptibility x vs. T/t, obtained by the high-T 
expansion. 

It is quite delicate to approach T = within the J = 
model. The ground state for a single hole is within the 




sector of maximal total spin S — Six 



181, however ex- 



cited states are numerous and close in energy, so that the 
transition between the regime of an incoherent hole prop- 
agation and the regime of the large FM polaron appears 
to happen at surprisingly low T* /t 0.1. This is consis- 
tent with the well established fact that the FM-polarized 
ground state is very sensitive to any change of parame- 
ters. The behavior of a single hole is simple only strictly 
at T = 0. It is expected that the QP at B = behaves ac- 
cording to the Nagaoka theorem, i.e. as a free hole with 
energy e = — 4 < in a filled band of spinless fermions. 



4 



Looking only in the sector S = 5'max, a finite field {B > 0) 
should increase the ground state energy according to the 
cyclotron frequency, i.e. 6e ^ eB/m* — Bt. 

In Fig. ^ we show both the absolute g.s. energy s and 
the lowest energy in the sector S = S'max as calculated 
within small systems with N = 16, 18 and 20. For the 
Nagaoka sector S = S'max the behavior is for a < n/A 
clearly of the cyclotron type, while for higher a there 
are some visible commensurability anomalies identical to 
the study of spinless fermions in a magnetic field [ p4[ . 
On the other hand, deviations of the absolute ground 
state from the naive result are much more pronounced. 
First, even the smallest B = Bq leads to the instability 
of the S = S'max state, and the actual spin of the g.s. is 
S < Smax- Nevertheless, 6e remains quite close to the 
cyclotron value. For higher B > Bq the ground state 
saturates quite abruptly to lowest spin S = 1/2 and g.s. 
energy e is much lower than in the S — Smax case. Note, 
however, that even an approximate validity of the simple 
cyclotron- frequency argument, Se oc \B\, implies a diver- 
gent susceptibility x(T-^O) oo, as found also from the 
high-T expansion. 



-2.8 



-3.2 



-3.6 



-4.0 



• 16 

.18 J=0 

♦ 20 

O 16 
□ 18 
O20 



4- ■■■ 



7t/4 



7l/2 



2ii/3 



FIG. 3. Ground state energy e (in units of t) vs. a, ob- 
tained via the diagonalization of small systems of different 
sizes and J = 0. Open dots correspond to the lowest-energy 
state within the sector S — 5'max. 



III. J > 0: ANTIFERROMAGNETIC SPIN 
BACKGROUND 

J > introduces quite a different QP behavior. At 
T = the spin background corresponds to an AFM with 
long range order. The ground state of the spin polaron is 
well understood, and corresponds within the t-J model 
to S = 1/2 and k = k*, with a weak dispersion along the 
AFM zone boundary. 

At least for weak magnetic fields it is expected that 
the QP description of the hole is still valid, resulting in 
a cyclotron motion with the linear-in-field dependence of 



the QP energy. Such behavior, is questioned by recent 
small-system diagonalisation results |^ . It should be re- 
minded that results obtained from the exact diagonali- 
sation can be misleading due to finite size effects, if the 
cyclotron radius is comparable to the linear size of the 
system. With p.b.c. even more pronounced effects com- 
pared to those already discussed for J = in Sec. II can 
be expected. For the AFM case they arise from larger 
ground state degeneracy, and from the near-degeneracy 
of the QP dispersion along the AFM zone boundary. 



A. Cumulant expansion method 

At T~0 we can consider a QP moving in an ordered 
AFM by performing the expansion starting in the limit 
t/J <^ 1. Here the standard cumulant-expansion (CE) 
procedure for the ground-state energy, as first considered 
by C. Bloch H, is followed. 

The implementation of the method in the zero-field 
case with a single hole in the t-J model has been already 
given in some detail elsewhere |l^. There the diagonal 
part of the Heisenberg spin interaction in Eq. (1) has 
been taken as the unperturbed Hamiltonian Hq and the 
ground-state energy £o(k), expressed relative to the un- 
doped AFM ground state, has been obtained as a double 
power series in u — t/J and 7 = J±/ J 



(k)^."( 



(17) 



where the k-dependence of the expansion coefficients 
merely refiects the underlying translational invariance of 
the problem. 

A finite magnetic field breaks the translational invari- 
ance of the system, so the above result, obtained from 
the nondegenerate theory is not valid. Indeed, the un- 
perturbed ground state becomes highly degenerate, al- 
lowing for any hole position in an otherwise perfectly 
Neel-ordered spin state, thus the degenerate perturba- 
tion theory must be used. The version of the secular de- 
terminant as given in is used, which up to the fourth 
order in perturbation H' reads 



= < ^h'9Rh'\j > + < i\H'^H'^H'^H'\] > 

m 



ji Qo Qo j^i Qo 
vo vo m 



< i\H'^H'\£ >< £\H'^H'\j > 



VO 



(18) 



The (bra)kets denote the set of degenerate states - po- 
sitions of the hole - with the unperturbed energy Eq, 
whereas Pq and Qq are the projectors onto this subset 
of states and its complement, respectively. H' is taken 
to be the sum of the hopping part and the transverse 
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spin part of H in (1), while 770 stands for the energy 
denominator Eq — Hq. A contribution to the matrix el- 
ement M{i,j) —< i\- ■ - {j > resulting from a particular 
process involving products of operators as in ( |l8| ) is 
calculated from the associated graph in any order 

of perturbation expansion. However, as the hole moves 
from some initial position i to a final but equivalent posi- 
tion j under the action of a particular string of operators 
Op along some path, it acquires a definite phase factor 
Oij, depending on the path itself. Two such paths are 
depicted in Fig. ^. 

In the chosen Landau gauge the phase associated with 
any link in the y-direction is given by the x-coordinate of 
that link. Thus, the total phase 9ij acquired by the hole 
along the path C, is 



iiC) 



A ■ ds ^ i^ijy -iy)a + (l)ij{C). (19) 



Here (j)ij{C) is the phase relative to the initial point at i. 
Care should be taken of the proper orientation in which 
the link is traversed. Thus, referring to Fig. ^, phase 
along Cp, is (j) — —2a, whereas along the path Af_i (j) = 4a. 
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FIG. 4. The phase increments, defined on the (oriented) 
links between two neighboring lattice points, in the Landau 
gauge A = B{0, x, 0) along two difi'erent paths and A^i 
connecting points i and j. 

We have generated all the paths in order u^j™ with 
(n,TO) up to (2,4), (4,4), (6,3) and (8,2). Each con- 
tribution to the perturbation series in order (n, m) is 
then given by the magnitude Lo{i,j) of the matrix ele- 



to a definite product of operators of order r = n + m 
along the path from i ^ j [1^. 

The secular equation (|l8) becomes a difference equa- 
tion for the on-site amplitudes fi on a rectangular grid 
with N — Lx X Ly sites 

e/, = ^Af(i,j)/„ 

j 

(20) 

where i and j run only over one sublattice, M{i,j) is a 
sum of contributions along different paths C^, and the 
energy e = e(a) is again measured with respect to the 
g.s. of the undoped AFM state. Referring again to Fig. ^ 
the paths and Af_i would then first appear in order 
u^-f^ and w'*7^, respectively. 

In the chosen gauge Eq. (|^), the system is translation- 
ally invariant along the y-direction. Thus, the ansatz 
fi — 9ix exp(ifcj^ij,) reduces the above equation to a dif- 
ference equation in one dimension, where £ — ix, 



1 < < 



Q — ky + a£, 



(21) 



and the summation over r is restricted to run over values 
for which £' — £ + t is even, whereas ky £ [0, 2tt]. Note 
also that lo and do not depend on the initial point £, 
but only on the path C. Taking a — 2Trp/Lx, one can 
impose the p.b.c. also in the x direction. 

The eigenvalue problem of Eq.(|2l|) is solved numeri- 
cally for Lx ^ where ^ = 13 is the smallest linear size 
of the region visited by the hole, within the order of per- 
turbation series here considered. In Fig. |^ we plot the g.s. 
energy e as a function of a at J/t = 2 and the isotropic 
exchange 7 = 1, evaluated for Lx — 128. A linear-in-field 
dependence of 6e{a) — e(a) — e(0) in Fig. || is evident for 
small a, implying that the hole may still be described 
as a QP exhibiting cyclotron motion. However, after the 
initial rise an almost monotonic decrease is observed and 
the minimum of s{a) is achieved for a = tt. The figure 
includes also the respective data from exact diagonaliza- 
tion for N — 20. Although both sets of data do not agree 
in detail, which can be partially attributed to the pertur- 
bational character of the CE results, the overall behavior 
is remarkably similar, including the initial rise in Se{a) 
for = 20 and the cusp-like behavior close to a = tt/2. 
This holds to a lesser degree for other commensurability 
points, e.g. for a = 7r/4 and 37r/4. 



ment {i\0^^^\j) and the phase 9^(i,j), where 0)1' refers 
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FIG. 5. Ground state energy difference 5e{a) = e{a) ~e(0) 
(in units of t), as obtained for t/J — 0.5 and 7 = 1 from the 
cumulant expansion method. For comparison small-system 
data for TV = 20 are also included (♦) 

The experimentally relevant value of J/t is ^ 0.3. 
Technically, it is possible to perform a Pade-like extrap- 
olation of s(a) to larger values of w, in analogy to the 
B — case . Still we do not attempt to perform such 
an extrapolation here. Relying on previous experience 
[0 that no crossover in the QP behavior sets in down to 
J/t <C 1 (where Nagaoka regime takes over) we believe 
that the qualitative conclusions on QP behavior for finite 
B remain valid in the physically relevant regime J < t as 
well. 



B. Small-system results 

Using the numerical finite temperature technique de- 
scribed in Sec. IIB and applied to the J — case the 
results for J > are also obtained. In the following 
we choose J/t = 0.4, which is close to the situation in 
cuprates, and consider the contribution of the single hole 
energy e = E{Nfi—l) — E{Nfi—0) to the internal en- 
ergy at T = and T > ((e)). First note that for 
T max(i, J), the leading order of the high-T expan- 
sion is independent of J and in this limit results for J = 
and J > match. 

In Fig. ^ the intermediate temperature hole energies 
(e(a)) are evaluated at two different T within different 
system sizes. Comparing results with those for J = 
in Fig. 0, several conclusions can be reached. At higher 
T/t = 1 energies {s{a)) for J = and J/t — 0.4 are 
qualitatively similar. As a function of a both cases cor- 
respond approximately to the simple cos a variation. It 
is however evident that finite J considerably reduces (by 
a factor ^ 4) the total energy span. 

Relative to the J — case it is also clear that finite 
size effects become more pronounced when approaching 
the low-T regime. This is related to the near-degeneracy 



of the quasiparticle dispersion eo(k) at the bottom of 
the QP band, which is poorly reproduced in small sys- 
tems. Nevertheless, numerical results establish quite con- 
sistently the crossover at T ^ T* (x J . 

For T < T* energy (e(a)) shows an overall opposite, 
i.e. a paramagnetic-like variation with a, as found in |^ 
and reproduced above within the CE with respect to t/J. 
The results in this temperature region are more size de- 
pendent. At T/t = 0.2 in Fig. | we can find a maximum 
for a > only in systems with N = 18, 20 sites, while the 
iV = 16 case shows a different behavior. The deviation 
within the latter system is due to additional symmetry 
of the 4x4 system (hypercube). 



A 
V 



-1.7 
-1.8 - 

A 

u 

^ -1.9 
-2.0 

-1.83 
-1.84 
-1.85 
-1.86 



1 1 1 

-0 


1 6 - 




• 18 




♦ 20 






■♦ ♦ ° 




♦ ° 








♦ 




• 




♦ 






•♦° 




•♦ 


T/t=0.4 


•♦ 




• • 


H 1 h- 

T/t=1 


1 h 






° • §" 














• ♦ 


. * 




• 




* 




• 












1 1 1 


J/t=0.4 



71/2 

a 



FIG. 6. Single hole energy (e) vs. a for J/t — 0.4 at two 
different T/t, as calculated on a system of N = 16, 18, 20 sites. 

Quite similar results are obtained for the ground state 
energy e, shown in Fig. 0. Again, N — 18 and iV = 20 
systems yield very similar e(a), while = 16 results de- 
viate especially for 'commensurate' values a = 0,7r/2,7r. 
The qualitative trend of e{a) agrees well with the CE 
results on Fig. 5. The similarity is in the existence of 
(rather shallow) maximum at low a for largest N, indi- 
cating a QP behavior with a cyclotron motion in weak B 
and also a pronounced reduction of e for a > 7r/2. The 
only difference is in the absence of commensurability dips 
in Fig. 8, which we attribute to the very particular sys- 
tem shapes used in the calculation, and to the almost 



7 



degenerate dispersion along the magnetic Brillouin zone 
boundary. 
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FIG. 7. Ground state energy e{a) for J/t 
lated on systems of A*' = 16, 18, 20 sites. 



0.4, as calcu- 



Finally, let us present in Fig. ^ results for the orbital 
susceptibility x{T), Eq.(^. For J = the calculation 
is performed via high-T expansion as already discussed 
in Sec. 2. Since we rely only on discrete values of a, the 
analysis of finite-system data for J/t — 0.4 is on the other 
hand less rehable. 
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FIG. 8. Susceptibility x vs. T/t for J/t — 0.4, as extracted 
from small-system data A'' = 18, 20. 

Moreover, the variation of F{a) is quite subtle at 
smallest amin, as evidenced also from e in Fig. 0. In 
Fig. _g| the susceptibility x is presented as obtained from 
with parabolic fit for F{a) using only a = and 
ctmin- This procedure could be questionable for 
T*, but qualitative behavior is still quite instruc- 
Relative to the J = case, the diamagnetism is 
suppressed by J > at higher T ^ J. In an intermedi- 
ate regime T ^ J, susceptibility x appears even to change 




sign, i.e. becomes paramagnetic. Only at low T < T* a 
pronounced strongly T-dependent diamagnetic response 
is again observed, consistent with the QP cyclotron mo- 
tion at T = 0. It is an interesting observation that an 
exact solution of the problem on a single plaquette (at 
T > 0) subjected to an effective staggered AFM field, 
reproduces qualitative features of Fig. ^. 



IV. CONCLUSIONS 

Our study shows that the calculation of the effects of 
finite magnetic field, coupled to orbital motion of elec- 
trons, becomes delicate in models of correlated electrons, 
e.g. within the t-J model considered in this work. Re- 
sults for magnetic observables such as the diamagnetic 
susceptibility appear to be strongly influenced by finite- 
size effects, which are hard to overcome in available sys- 
tem sizes. Both for J = and J > some of deviations 
appear at 'commensurate' values of a within the given 
system geometry, and are particularly large for iV = 16. 
We have shown that it is possible to understand such 
finite-size effects within the high-T expansion or within 
the t/J perturbation expansion as a contribution of addi- 
tional graphs due to p.b.c. Still it is impossible to elim- 
inate them systematically in most interesting physical 
regimes. These effects lead to a nonmonotonous varia- 
tion of observables, e.g. {s{a)), which in turn leads to an 
enhanced uncertainty in x(T). 

The J = case seems both easier to study and to 
understand. High-T expansion and small systems show 
a continuous transition from the high-T regime of in- 
coherent hopping to the Nagaoka state at T = 0, with 
a monotonous increase of the diamagnetic x- At T > 
where the HTE is reliable the variation of the energy with 
field {s{a)) is quite close to a simple cos a form. Never- 
theless the asymptotic behavior at low T is not simple to 
establish, since the nature of low lying states (above the 
Nagaoka state) is complicated. 

The behavior of the AFM J > polaron is more in- 
volved. While at T ^ J the exchange scale J is not im- 
portant and results qualitatively follow those for J = 0, 
new physics emerges for T < J. A nearly flat (^(q;)) at 
intermediate regime T '--^ J is quite remarkable, and leads 
to a vanishing diamagnetic x (or even change of its sign) . 
It seems, that here J > diminishes and even destroys 
emerging coherence of QP. Only at lower T < J the co- 
herence is established and the known dynamical picture 
of a coherent AFM polaron is dominating the behavior 
in lowest fields a. Reliable results are however difficult 
to obtain even for T = 0, since in small systems their 
variation with a is very sensitive to the system shape, 
boundary conditions etc. due to the very anisotropic and 
degenerate QP dispersion. 

The CE results are very instructive, but it is not 
straightforward to make them quantitative for J < t. 
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For J > t, where the CE series converges quite rapidly, 
there is a clear evidence in weak fields that the QP is 
exhibiting cyclotron motion. The commensurability ef- 
fects are quite pronounced and agree with finite-cluster 
data for large J/t > 1, see Fig. |5[ On the other hand, 
these effects are not evident in small-system results for 
J/t < 1. This non- agreement is attributed in part to 
specific lattice shapes and p.b.c. The flatness of the dis- 
persion seems to be the main reason for the structureless 
character of small-system data in the region of small a 
for J/t<l. 

Here we would like to point to the close similarity with 
the challenging theoretical problem of the Hall effect, 
which is difhcult to approach even for very low doping, 
e.g. for a single hole |^8|. Similarly to the orbital suscep- 
tibility the Hall effect emerges due to coupling to orbital 
currents. The Hall constant is given by Rh — f^xy/Ba^^, 
where the off-diagonal conductivity axy can be related to 
the orbital susceptibility x as H] 



J ^ B— — 

dch d^i ' 



(22) 



where /i is the chemical potential. In our case of very 
low hole doping, i.e. in the semiconductor-like regime, 
the susceptibility x scales linearly with Ch] moreover 
dch/d^ = —(3ch, so that axy oc x- The high-T expansion 
of the Hall constant R*jj{T) (the high-frequency value) 
is analogous to that of x(T) On the other hand, 
crossing the scale T ^ J remains the challenge, whereby 
it seems that at this intermediate T the hole-like R*jj is 
even reduced with respect to its high-T value H,^. Ex- 
periments 0] indicate, that Rh recovers for T < T*, 
varying strongly with T and approaches the well known 
quasiclassical result for T ^ |2^ . Note that our results 
for xC^); Fig. ^, indicate just on such behavior. 

Let us finally comment on the magnitude of the dia- 
magnetic susceptibility. Since we are evaluating the case 
of a single hole, at low doping c/j ^ 1 the observable dia- 
magnetic contribution to the susceptibility (per unit cell) 
should be X = Cc/jXO) where ( is dimensionless value, pre- 
sented in Figs. ^ and || It is convenient to compare these 
values to the spin susceptibility (per unit cell) of the pla- 
nar undoped AFM for T < T* , where Xs ~ i-OfiQ^is/J 
1^,^. Setting rut = /2taQ the ratio can be expressed 
as 



X_ 

Xs 



me 
nit 



(23) 



where K 2.8 is a numerical constant. Taking the stan- 
dard values for cuprates t = 0.4 eV, J/t = 0.3 and oq = 
0.38 nm, the above relation reduces to x/Xs ~ 1-9 C.. 

Another relation can be obtained for the Pauli suscep- 
tibility of a half-filled band of free tight-binding electrons, 
where a constant (average) density of states is assumed 
for simplicity. This gives a similar relation 



^=/^'CcJ^) ^9.3c,C, (24) 

XP ymtj 

with K' ^ 4.0. To estimate the actual value of x/Xs in 
Eq. (27) we take Ch ^ 0.15, e.g. as in the 'optimal' dop- 
ing regime, and C —0.1 in the region T < T*. This 
gives x/Xs ~ —0.03, which is of the same order of mag- 
nitude as the experimentally measured value [Q. Note, 
however, that below the crossover temperature T* C be- 
comes strongly temperature dependent, as opposed to 
the usual T-independent Landau-type diamagnetism in 
Fermi liquids. Since it is difficult to distinguish different 
contributions to the actual susceptibility in experiments, 
it remains to be seen whether such T-dependent x really 
appears in cuprates and analogous systems. 
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